Random number generation method selecting system, random number generation method selecting method, and random number generation method selecting program

ABSTRACT

A generation means  11  generates a uniform random number between 0 and a first probability, which is a probability of a stochastic variable becoming a value within a predetermined interval in a positive range in the first discrete distribution. When a uniform random number less than or equal to a second probability is generated, the second probability being a probability of the stochastic variable becoming a value within a predetermined interval in a second discrete distribution, which is a discrete Gaussian distribution on a one-dimensional lattice the center of which is the origin, the selection means  12  selects, as a random number generation method, an accumulation method in which a functional value defining the second discrete distribution is used. When a uniform random number greater than the second probability is generated, the selection means  12  selects a rejection sampling method as the random number generation method.

This application is a National Stage Entry of PCT/JP2017/028584 filed onAug. 7, 2017, the contents of all of which are incorporated herein byreference, in their entirety.

TECHNICAL FIELD

The present invention relates to a random number generation methodselecting system, a random number generation method selecting method,and a random number generation method selecting program, andparticularly relates to a random number generation method selectingsystem, a random number generation method selecting method, and a randomnumber generation method selecting program that are used inlattice-based cryptography and signature and generates a random numberaccording to discrete Gaussian distribution the center of which is notan origin.

BACKGROUND ART

First, discrete Gaussian distribution will be defined. A functiondefined by a real number s∈R (R is a symbol representing a set of allreal numbers) is defined as follows.

$\begin{matrix}{\left\lbrack \text{Expression~~1} \right\rbrack\mspace{490mu}} & \; \\{{\phi_{s}(x)}:={\frac{1}{s}{\exp\left( {{- \frac{\pi}{s^{2}}}x^{2}} \right)}}} & {{Formula}\mspace{14mu}(1)}\end{matrix}$

A distribution in which an integer value u∈Z (Z is a symbol representinga set of whole integers) is output with probability φ_(s)(u)/Σ^(∞)_(j=−∞)φ_(s)(j) is referred to as discrete Gaussian distribution with avariance value s. As described in Non Patent Literature (NPL) 1, arandom number generated in accordance with the above discrete Gaussiandistribution is used for cryptography using a lattice (hereinafterreferred to as lattice-based cryptography). Lattice-based cryptographyis expected to be used as post-quantum cryptography. Furthermore,lattice-based cryptography is a cryptographic system that has beenstudied as a cryptographic scheme with high computational efficiency andhigh functionality.

The discrete Gaussian distribution is probability distribution that canoutput all integer values. However, when the random number generated inaccordance with discrete Gaussian distribution are used forlattice-based cryptography, limiting the range of integer values outputby the discrete Gaussian distribution would often increase efficiency inthe generation of stochastic variables.

For example, limiting the range of integer values output by discreteGaussian distribution in such a way that it depends on a securityparameter n∈N (N is a symbol representing a set of all natural numbers)would increase the efficiency in the generation of stochastic variables.

That is, setting {k∈Z|−t·s≤k≤t·s} as an output range of the discreteGaussian distribution with t=ω(logn)^(1/2) would increase the efficiencyin the generation of stochastic variables. Here, ω is a Landau symbol.It is generally known that limiting the range of integers output bydiscrete Gaussian distribution as described above would not affect thesecurity of lattice-based cryptography.

When the normalization constant W is defined as W=Σ^(t·s)_(i=−t·s)φ_(s)(i), the discrete Gaussian distribution in which the rangeof output integers is limited as described above is expressed asΨ_(s)(x)=φ_(s)(x)/W. When Ψ_(s)(x) is used, the integer value u∈Z isoutput with probability Ψ_(s)(u).

Hereinafter, the “discrete Gaussian distribution” in this specificationwill refer to probability distribution in which an output integer rangeis {k∈Z|−t·s≤k≤t·s}, and an integer value u∈Z is output with probabilityΨ_(s)(u). Moreover, the function Ψ_(s)(x)=φ_(s)(x)/W will be referred toas a function that defines discrete Gaussian distribution.

Next, the center of the discrete Gaussian distribution will bedescribed. The discrete Gaussian distribution with the center c and thevariance value s is probability distribution that outputs an integervalue u with probability Ψ_(s)(u−c).

The above is a definition for discrete Gaussian distribution on aone-dimensional lattice. Next, a definition for discrete Gaussiandistribution on an n-dimensional lattice will be given.

A matrix in which vectors {b₁ ^(→), . . . , b_(n) ^(→)}∈R^(n) arearranged horizontally is denoted as B. An n-dimensional lattice A(B)using the matrix B is defined as follows.

$\begin{matrix}{\left\lbrack \text{Expression~~2} \right\rbrack\mspace{490mu}} & \; \\{{\Lambda(B)}:=\left\{ {\sum\limits_{i\; \in {\lbrack n\rbrack}}{c_{i} \cdot {\overset{\_}{b}}_{i}}} \middle| {c_{i} \in \; Z} \right\}} & {{Formula}\mspace{14mu}(2)}\end{matrix}$

The Gaussian function on R^(n) whose center is c^(→) is defined asfollows using the parameter s.

$\begin{matrix}{\left\lbrack \text{Expression~~3} \right\rbrack\mspace{490mu}} & \; \\{{{{for}\mspace{14mu}{all}\mspace{14mu}\overset{\rightarrow}{x}} \in R^{n}},{{\rho_{s,\overset{\rightarrow}{c}}\left( \overset{\rightarrow}{x} \right)} = {\exp\left( {{- \pi}\frac{{{\overset{\rightarrow}{x} - \overset{\rightarrow}{c}}}^{2}}{s^{2}}} \right)}}} & {{Formula}\mspace{14mu}(3)}\end{matrix}$

The discrete Gaussian distribution on the n-dimensional lattice A isdefined as follows using the Gaussian function expressed by the aboveFormula (3).

$\begin{matrix}{\left\lbrack \text{Expression~~4} \right\rbrack\mspace{490mu}} & \; \\{{{for}\mspace{14mu}{all}\mspace{14mu}{\overset{\rightarrow}{x} \in \Lambda}},{{D_{\Lambda,s,\overset{\rightarrow}{c}}\left( \overset{\rightarrow}{x} \right)} = \frac{\rho_{s,\overset{\rightarrow}{c}}\left( \overset{\rightarrow}{x} \right)}{\rho_{s,\overset{\rightarrow}{c}}(\Lambda)}}} & {{Formula}\mspace{14mu}(4)}\end{matrix}$

That is, the discrete Gaussian distribution on the n-dimensional latticeis probability distribution in which the stochastic variable follows theFormula (4). Hereinafter, for simplicity, the discrete Gaussiandistribution on the one-dimensional lattice and the discrete Gaussiandistribution on the n-dimensional lattice will also be referred to asone-dimensional discrete Gaussian distribution and an n-dimensionaldiscrete Gaussian distribution, respectively.

The above is the definition of one-dimensional discrete Gaussiandistribution and the definition of a multidimensional (n-dimensional)discrete Gaussian distribution. Next, a sampling method for generating arandom number according to each of discrete Gaussian distributionpatterns will be described.

Typically, sampling methods as a method for generating a random numberaccording to one-dimensional discrete Gaussian distribution include twomethods, an accumulation method and a rejection sampling method. Here, afunction that defines one-dimensional discrete Gaussian distribution isφ(x), and an output range of the one-dimensional discrete Gaussiandistribution is {k∈Z|−t·s≤k≤t·s}.

The above two sampling methods will be described separately for the casewhere a center of the one-dimensional discrete Gaussian distribution isan origin and the case where a center of the one-dimensional discreteGaussian distribution is not the origin.

First, a process of generating, using an accumulation method, a randomnumber according to one-dimensional discrete Gaussian distribution thecenter of which is the origin will be described on the basis of thedescription of NPL 2. FIG. 9 is a block diagram showing an exemplaryconfiguration of a conventional random number generation system thatgenerates a random number according to one-dimensional discrete Gaussiandistribution.

As shown in FIG. 9, a random number generation system 910 that generatesa random number according to one-dimensional discrete Gaussiandistribution the center of which is the origin includes a storage device911, a searcher 912, and a uniform random number generator 913.

The following describes operation of the conventional random numbergeneration system 910 including components as shown in FIG. 9 thatgenerates a random number according to one-dimensional discrete Gaussiandistribution the center of which is the origin. First, values of φ(0)/2,φ(1), φ(2), . . . , φ(t·s) are computed in advance and stored in thestorage device 911. Individual vertical lines in the storage device 911shown in FIG. 9 represent individual values stored, such as φ(0)/2 orφ(1).

Next, the uniform random number generator 913 outputs a real valuex∈[0,1]. The output x∈R is input to the searcher 912. The searcher 912to which x has been input performs binary search for z∈Z that satisfiesφ(z−1)≤x<φ(z) from among the values stored in the storage device 911.

Next, the searcher 912 uniformly selects a code sign=±. Next, thesearcher 912 outputs sign·z∈Z as a random number according toone-dimensional discrete Gaussian distribution. The above is a method,using the accumulation method, for generating a random number accordingto one-dimensional discrete Gaussian distribution with the origin at thecenter.

Note that the method for generating the random number according to theone-dimensional discrete Gaussian distribution the center of which isnot the origin by the accumulation method, for example, the center beingat a, can be a method to replace φ(x) with φ(x−a) in the method forgenerating the random number according to the one-dimensional discreteGaussian distribution the center of which is the origin.

In the accumulation method, the number of pieces of data stored in thestorage device is proportional to t·s. The data stored in the storagedevice when using the one-dimensional discrete Gaussian distribution thecenter of which is the origin would be φ(0)/2, φ(1), φ(2), . . . ,φ(t·s). In contrast, the data when using one-dimensional discreteGaussian distribution the center of which is a, rather than the origin,would be φ(0−a)/2, φ(1−a), φ(2−a), . . . , φ(t·s−a).

That is, the accumulation method has a problem that using discreteGaussian distribution having a large variance value s would increase theamount of memory required to store the functional values. In theexisting lattice-based cryptography, the variance value s takes arelatively large value.

As a specific amount of memory required for storing the functionalvalues, for example, the values shown in FIG. 10 are described in NPL 4.FIG. 10 is an explanatory diagram showing an example of the amount ofmemory consumed when the accumulation method is used. As shown in FIG.10, the larger the variance value, the larger the amount of memoryconsumed.

The description of “q-type/1 signature” shown in the “center” field ofFIG. 10 indicates that the center of q=2^(K) will be prepared for onesignature generation. Similarly, “2-type/1 signature” indicates thatonly two types of centers will be prepared for one signature generation.

In addition, “Usage” shown in FIG. 10 represents an application in whichGaussian distribution is used by individual methods. For example, in themethod “LWE-plane”, Gaussian distribution is used for the purpose of“key generation”. Furthermore, in the method “GPV” or the like, Gaussiandistribution is used for the purpose of “signature generation”. WhenGaussian distribution is used in the “signature generation” application,the Gaussian distribution is used many times.

In general, lattice-based cryptography such as RSA requires a smallamount of computation. Therefore, the lattice-based cryptography isexpected to be used in a device having a small computing resource and asmall storage capacity, such as a sensor device or a mobile phone.

However, when the accumulation method is used for generation of a randomnumber according to the discrete Gaussian distribution, which is alattice-based cryptography subroutine, a large amount of storagecapacity would be used for storing data used in the accumulation method.That is, there would be a problem of difficulty using the lattice-basedcryptography in a device having a small storage capacity.

Next, the rejection sampling method will be described. The rejectionsampling method is a method used for generation of a random numberaccording to any discrete probability distribution, not limited to thediscrete Gaussian distribution.

First, a rejection sampling method for generating a random numberaccording to discrete probability distribution for a general stochasticvariable X will be described on the bases of description of NPL 5.Thereafter, a rejection sampling method for generating a random numberaccording to one-dimensional discrete Gaussian distribution will bedescribed.

In order to generate a random number according to discrete probabilitydistribution p(X=x_(i)) using the rejection sampling method, a functiont(x) that can be efficiently computed is prepared, from among thefunctions t(x) that satisfy t(x_(i))≥p(x_(i)) for all x_(i) values.

Next, a function r(x) in which t(x) is normalized is set asr(x)=t(x)/Σt(x_(i)). Next, the following procedure is executed togenerate a random number according to the discrete probabilitydistribution p(X=x_(i)) (probability distribution function) for thestochastic variable X by using the rejection sampling method.

(Step 1) Generating a random number Y according to the probabilitydistribution function r(x).

(Step 2) Generating a uniform random number U in an interval [0,1]independently of Y.

(Step 3) When U≤p(Y)/t(Y) is satisfied, the random number X is set toX=Y. When U≤p(Y)/t(Y) is not satisfied, the process of (Step 1) will beperformed again.

Note that 1/ε, g, and f in Algorithm A₁ described in NPL 5 are convertedto 1, t, and p, respectively, in the above procedure.

Next, a rejection sampling method for generating a random numberaccording to one-dimensional discrete Gaussian distribution will bedescribed on the basis of description of NPL 3. Note that the methoddescribed in NPL 3 is a method in which the above-described rejectionsampling method is applied to a case where the function t(x) isidentically 1.

FIG. 11 shows an exemplary configuration of a device for implementingthe method described in NPL 3. FIG. 11 is a block diagram showinganother exemplary configuration of a conventional random numbergeneration system that generates a random number according toone-dimensional discrete Gaussian distribution.

As shown in FIG. 11, a random number generation system 920 thatgenerates a random number according to one-dimensional discrete Gaussiandistribution the center of which is the origin includes a rejectiondeterminer 921, a uniform random number generator 922, and an outputdevice 923.

The following describes operation of the conventional random numbergeneration system 920 including components as shown in FIG. 11 thatgenerates a random number according to one-dimensional discrete Gaussiandistribution the center of which is the origin. First, the uniformrandom number generator 922 generates a uniform random number u₁∈Z in arange of {k∈Z|−t·s≤k≤t·s}.

Next, the uniform random number generator 922 generates a real-valuedrandom number u₂∈R within a range of [0, φ(0)]. The uniform randomnumber generator 922 inputs the generated u₁∈Z and u₂∈R to the rejectiondeterminer 921.

Next, rejection determiner 921 compares φ(u₁) with u₂∈R. When thecomparison result is u₂≤φ(u₁), the rejection determiner 921 inputs u₁∈Zto the output device 923. The output device 923 outputs u₁∈Z as a randomnumber according to one-dimensional discrete Gaussian distribution thecenter of which is the origin.

When the comparison result is u₂>φ(u₁), the random number generationsystem 920 returns to the first step and executes the same operationagain. The above is a method, using the rejection sampling method, forgenerating a random number according to one-dimensional discreteGaussian distribution with the origin at the center.

Note that the method for generating the random number according to theone-dimensional discrete Gaussian distribution the center of which isnot the origin by the rejection sampling method, for example, the centerbeing at a, can be a method to replace φ(x) with φ(x−a) in the methodfor generating the random number according to the one-dimensionaldiscrete Gaussian distribution the center of which is the origin.

In the rejection sampling method, when the condition in the processcorresponding to the above (Step 3) is not satisfied, the similarprocess will be repeatedly executed. That is, the rejection samplingmethod is required to recompute a function that defines the discreteGaussian distribution every time the generated uniform random number isrejected, leading to a problem of reduction in the computationefficiency.

In summary, the accumulation method has an advantage that thecomputation cost is low. On the other hand, the accumulation method hasa disadvantage that the memory cost is high. The rejection samplingmethod has an advantage that the memory cost is low. On the other hand,the rejection sampling method has a disadvantage that the computationcost is high.

Next, a method for generating a random number according to discreteGaussian distribution on a multidimensional lattice will be described onthe basis of the description of NPL 3. The following are those preparedbefore the explanation.

Each of Gram-Schmidt orthogonalized vectors a₁ ^(˜→), . . . , a_(n)^(˜→) for vectors a₁ ^(→), . . . , a_(n) ^(→) is set to a vectorcomputed as follows.

$\begin{matrix}{{{\overset{\rightarrow}{\overset{˜}{a}}}_{1} = {\overset{\rightarrow}{a}}_{1}}{{\overset{\rightarrow}{\overset{˜}{a}}}_{2} = {{\overset{\rightarrow}{a}}_{2} - {\frac{\left\langle {{\overset{\rightarrow}{a}}_{2},{\overset{\rightarrow}{\overset{˜}{a}}}_{1}} \right\rangle}{{{\overset{\rightarrow}{\overset{˜}{a}}}_{1}}^{2}}{\overset{\rightarrow}{\overset{˜}{a}}}_{1}}}}{{\overset{\rightarrow}{\overset{˜}{a}}}_{3} = {{\overset{\rightarrow}{a}}_{3} - {\frac{\left\langle {{\overset{\rightarrow}{a}}_{3},{\overset{\rightarrow}{\overset{˜}{a}}}_{2}} \right\rangle}{{{\overset{\rightarrow}{\overset{˜}{a}}}_{2}}^{2}}{\overset{\rightarrow}{\overset{˜}{a}}}_{2}} - {\frac{\left\langle {{\overset{\rightarrow}{a}}_{3},{\overset{\rightarrow}{\overset{˜}{a}}}_{1}} \right\rangle}{{{\overset{\rightarrow}{\overset{˜}{a}}}_{1}}^{2}}{\overset{\rightarrow}{\overset{˜}{a}}}_{1}}}}\vdots{{\overset{\rightarrow}{\overset{˜}{a}}}_{n} = {{\overset{\rightarrow}{a}}_{n} - {\sum\limits_{i = 1}^{n - 1}{\frac{\left\langle {{\overset{\rightarrow}{a}}_{n},{\overset{\rightarrow}{\overset{˜}{a}}}_{i}} \right\rangle}{{{\overset{\rightarrow}{\overset{˜}{a}}}_{i}}^{2}}{\overset{\rightarrow}{\overset{˜}{a}}}_{i}}}}}} & \left\lbrack {{Expression}\mspace{14mu} 5} \right\rbrack\end{matrix}$

In this specification, symbols “−”, “→”, “˜”, etc., which are symbolsused in the text should be described just above the previous character.However, due to restrictions on the text notation, these symbols will bewritten immediately after the character as described above. In theformulas and drawings, these symbols are written in their originalpositions.

FIG. 12 is an explanatory diagram showing an example of a random numbergeneration algorithm according to discrete Gaussian distribution on amultidimensional lattice. In the loop of the fourth to ninth lines ofthe algorithm shown in FIG. 12, a random number z_(i) is generated byone process, and c_(i−1) is updated on the basis of the generated z_(i).In the next processing, a random number z_(i−1) is generated after thediscrete Gaussian distribution is updated on the basis of the updatedc_(i−1). That is, a random number z_(n), . . . , z₁ are generatedsequentially.

In the process on the seventh line of the algorithm shown in FIG. 12,D_(Z,c′i,σi) represents one-dimensional discrete Gaussian distributionhaving a center c′_(i) and a variance value σ₁. Since c_(i−1) is updatedin each of processes that constitute the loop, generation of a randomnumber according to the one-dimensional discrete Gaussian distributionwith different centers is also required every time.

CITATION LIST Non Patent Literature

NPL 1: Regev, “On lattices, learning with errors, random linear codes,and cryptography,” STOC 2005, ACM, 2005, pages 84-93.

NPL 2: Chris Peikert, “An efficient and parallel Gaussian Sampler forlattices,” CRYPTO, 2010, pages 80-97.

NPL 3: Craig Gentry, Chris Peikert, and Vinod Vaikuntanathan, “How toUse a Short Basis: Trapdoors for Hard Lattices and New CryptographicConstructions,” STOC, 2008, pages 197-206.

NPL 4: DWARAKANATH, N. C, GALBRAITH, S. D, “Sampling From DiscreteGaussians for Lattice-Based Cryptography on a Constrained Device,” Appl.Algebra Engrg. Comm. Comput. 25, 2014, pages 159-180.

NPL 5: George Casella, Christian P. Robert, and Martin T. Wells,“Generalized Accept-Reject sampling schemes,” A Festschrift for HermanRubin Institute of Mathematical Statistics Lecture Notes—MonographSeries Vol. 45, 2004, pages 342-347.

SUMMARY OF INVENTION Technical Problem

In the case of generating a random number according to discrete Gaussiandistribution on a multidimensional (n-dimensional) lattice, it isrequired to perform n times of generation of a random number accordingto discrete Gaussian distribution on a one-dimensional lattice with acenter being not necessarily at the origin. The above accumulationmethod and rejection sampling method are used as a method of generatinga random number according to discrete Gaussian distribution on aone-dimensional lattice with a center being not necessarily at theorigin.

A problem when generating a random number according to the discreteGaussian distribution on the multidimensional lattice using theaccumulation method and a problem when generating a random numberaccording to the discrete Gaussian distribution on the multidimensionallattice using the rejection sampling method will be described in thisorder.

In a case where a random number is generated in accordance with thediscrete Gaussian distribution on the multidimensional lattice by usingthe accumulation method, a random number according to theone-dimensional discrete Gaussian distribution is generated by thenumber of dimensions of the lattice as an output target. In addition,the accumulation method is required to store a new numerical value inthe storage device for every different center of the one-dimensionaldiscrete Gaussian distribution used for generating a random number.

Therefore, when there is a substantial difference in each of the centersof a plurality of one-dimensional discrete Gaussian distributions usedfor generating a random number, a large amount of memory would beconsumed. That is, the accumulation method used for generating a randomnumber according to discrete Gaussian distribution on a multidimensionallattice would not be an efficient sampling method in terms of the amountof memory.

A problem when generating a random number according to the discreteGaussian distribution on a multidimensional lattice using the rejectionsampling method is that the generation speed of a random numberaccording to the one-dimensional discrete Gaussian distribution is low.The reason is that, every time the generated uniform random number isrejected, computation of a function that defines discrete Gaussiandistribution is required, lowering the overall computation efficiency.

Object of the Invention

In view of the above, the present invention provides a random numbergeneration method selecting system, a random number generation methodselecting method, and a random number generation method selectingprogram, for solution of the problems, that can further reduce thememory cost and the computation cost for generating a random numberaccording to discrete Gaussian distribution on a multidimensionallattice.

Solution to Problem

A random number generation method selecting system according to thepresent invention is a random number generation method selecting systemthat generates a random number according to a first discretedistribution, which is a discrete Gaussian distribution on aone-dimensional lattice the center of which is a positive value, therandom number generation method selecting system including: a generationmeans which generates a uniform random number between 0 and a firstprobability, which is a probability of a stochastic variable becoming avalue within a predetermined interval in a positive range in the firstdiscrete distribution; and a selection means, when a uniform randomnumber less than or equal to a second probability is generated, thesecond probability being a probability of the stochastic variablebecoming a value within a predetermined interval in a second discretedistribution, which is a discrete Gaussian distribution on aone-dimensional lattice the center of which is the origin, selects, as arandom number generation method, an accumulation method in which afunctional value defining the second discrete distribution is used, whena uniform random number greater than the second probability isgenerated, selects a rejection sampling method as the random numbergeneration method.

A random number generation method selecting method according to thepresent invention is a random number generation method selecting methodto be executed in a random number generation method selecting systemthat generates a random number according to a first discretedistribution, which is a discrete Gaussian distribution on aone-dimensional lattice the center of which is a positive value, therandom number generation method selecting method including: generating auniform random number between 0 and a first probability, which is aprobability of a stochastic variable becoming a value within apredetermined interval in a positive range in the first discretedistribution; when a uniform random number less than or equal to asecond probability is generated, the second probability being aprobability of the stochastic variable becoming a value within apredetermined interval in a second discrete distribution, which is adiscrete Gaussian distribution on a one-dimensional lattice the centerof which is the origin, selecting, as a random number generation method,an accumulation method in which a functional value defining the seconddiscrete distribution is used; and when a uniform random number greaterthan the second probability is generated, selecting a rejection samplingmethod as the random number generation method.

A random number generation method selecting program according to thepresent invention is a random number generation method selecting programthat is executed on a computer that generates a random number accordingto a first discrete distribution, which is a discrete Gaussiandistribution on a one-dimensional lattice the center of which is apositive value, the random number generation method selecting programcausing the computer to execute: a generation process that generates auniform random number between 0 and a first probability, which is aprobability of a stochastic variable becoming a value within apredetermined interval in a positive range in the first discretedistribution; a first selection process that, when a uniform randomnumber less than or equal to a second probability is generated, thesecond probability being a probability of the stochastic variablebecoming a value within a predetermined interval in a second discretedistribution, which is a discrete Gaussian distribution on aone-dimensional lattice the center of which is the origin, selects, as arandom number generation method, an accumulation method in which afunctional value defining the second discrete distribution is used; anda second selection process that, when a uniform random number greaterthan the second probability is generated, selects a rejection samplingmethod as the random number generation method.

Advantageous Effects of Invention

According to the present invention, it is possible to further reduce thememory cost and the computation cost for generating a random numberaccording to discrete Gaussian distribution on a multidimensionallattice.

BRIEF DESCRIPTION OF DRAWINGS

FIG. 1 is an explanatory diagram showing one-dimensional discreteGaussian distribution the center of which is an origin andone-dimensional discrete Gaussian distribution the center of which isnot the origin.

FIG. 2 is an explanatory diagram showing a region to which anaccumulation method is applied and a region to which a rejectionsampling method is applied.

FIG. 3 is a block diagram showing an exemplary configuration of a randomnumber generation system according to a first exemplary embodiment ofthe present invention.

FIG. 4 is a block diagram showing an exemplary configuration of ageneration method selection device 110 according to the first exemplaryembodiment.

FIG. 5 is a block diagram showing an exemplary configuration of arejection sampling device 120 according to the first exemplaryembodiment.

FIG. 6 is a block diagram showing an exemplary configuration of anaccumulation method sampling device 130 according to the first exemplaryembodiment.

FIG. 7 is a flowchart showing an operation of a random number generationprocess performed by a random number generation system 100 of the firstexemplary embodiment.

FIG. 8 is a block diagram showing a summary of a random numbergeneration system according to the present invention.

FIG. 9 is a block diagram showing an exemplary configuration of aconventional random number generation system that generates a randomnumber according to one-dimensional discrete Gaussian distribution.

FIG. 10 is an explanatory diagram showing an example of the amount ofmemory consumed when the accumulation method is used.

FIG. 11 is a block diagram showing another exemplary configuration of aconventional random number generation system that generates a randomnumber according to one-dimensional discrete Gaussian distribution.

FIG. 12 is an explanatory diagram showing an example of a random numbergeneration algorithm according to discrete Gaussian distribution on amultidimensional lattice.

DESCRIPTION OF EMBODIMENTS First Exemplary Embodiment

Hereinafter, exemplary embodiments of the present invention will bedescribed with reference to the drawings. In the present exemplaryembodiment, it is assumed that one-dimensional discrete Gaussiandistribution the center of which is the origin is virtually convertedinto one-dimensional discrete Gaussian distribution the center of whichis not the origin by the following method.

FIG. 1 is an explanatory diagram showing one-dimensional discreteGaussian distribution the center of which is the origin andone-dimensional discrete Gaussian distribution the center of which isnot the origin. For simplification, FIG. 1 shows only a positive rangeof an output range of the one-dimensional discrete Gaussiandistribution. As shown in FIG. 1, since the one-dimensional discreteGaussian distribution is a line-symmetric distribution, it is sufficientas long as the positive range of the output range is considered.

Referring to FIG. 1, it is observed that most of the one-dimensionaldiscrete Gaussian distribution with a center at a rather than at theorigin is included in the one-dimensional discrete Gaussian distributionthe center of which is the origin. For example, probability φ(1) in afirst interval shown in FIG. 1 (for convenience, the length of thehorizontal axis is 1) covers most of probability φ(1−a).

Similarly, probability φ(2) in a second interval shown in FIG. 1 (forconvenience, the length of the horizontal axis is 1) covers most ofprobabilityφ(2−a). Hereinbelow, similar characteristics can be observedin all intervals.

That is, the random number in the first interval according to theone-dimensional discrete Gaussian distribution the center of which isnot the origin is considered to be a random number according to theone-dimensional discrete Gaussian distribution the center of which isthe origin with the probability of φ(1)/φ(1−a). The present exemplaryembodiment utilizes the above characteristic.

FIG. 2 is an explanatory diagram showing a region to which anaccumulation method is applied and a region to which a rejectionsampling method is applied. Probability P shown in FIG. 2 is a sum ofthe probability φ(1−a), . . . , probability φ (t·s−a). That is, theprobability P is probability that a random number is output inaccordance with one-dimensional discrete Gaussian distribution thecenter of which is not the origin.

In addition, as shown in FIG. 2, probability p_({1-dim}) is a sum ofprobability φ(1), . . . , probability φ(t·s). That is, the probabilityp_({1-dim}) is probability that a random number is output in accordancewith one-dimensional discrete Gaussian distribution the center of whichis the origin.

In a case where a uniform random number r taken in an interval [0, P] isp_({1-dim}) or less, the random number generation system according tothe present exemplary embodiment generates a random number according tothe one-dimensional discrete Gaussian distribution the center of whichis not the origin by the accumulation method using a value of thefunction φ(x) that defines one-dimensional discrete Gaussiandistribution the center of which is the origin.

Furthermore, as shown in FIG. 2, probability p_({rjc}) is a sum ofprobability {(φ (1−a)−φ(1)}, . . . , probability {φ(t·s−a)−φ(t·s)}. Thatis, the probability p_({rjc}) is a difference between the probability Pthat a random number is output in accordance with one-dimensionaldiscrete Gaussian distribution the center of which is not the origin andthe probability p_({1-dim}) that a random number is output in accordancewith one-dimensional discrete Gaussian distribution the center of whichis the origin.

In a case where the uniform random number r taken in the interval [0, P]is larger than p_({1-dim}), the random number generation system of thepresent exemplary embodiment generates, by the rejection samplingmethod, the random number according to the one-dimensional discreteGaussian distribution the center of which is not the origin.

Since the probability p_({1-dim}) is sufficiently larger than theprobability p_({rjc}), the probability that a random number will begenerated by the rejection sampling method is low. That is, since thechance of generating a random number by the rejection sampling method isreduced, the overall computation cost for generating a random number isreduced.

In addition, since the accumulation method always uses the value of thefunction φ(x) that defines the one-dimensional discrete Gaussiandistribution the center of which is the origin, it is sufficient as longas the value of the function φ(x) is stored in the storage device. Thatis, the overall memory cost for generating a random number will also bereduced.

The above has described the method of generating a random numberaccording to one-dimensional discrete Gaussian distribution having thecenter shifted from the origin to the right. However, a random numberaccording to one-dimensional discrete Gaussian distribution having thecenter shifted from the origin to the left are generated in a similarmanner. The random number generation system according to the presentexemplary embodiment virtually converts the center of the discreteGaussian distribution on the one-dimensional lattice and therebyefficiently generates a random number according to the discrete Gaussiandistribution on the multidimensional lattice.

Description of Configuration

FIG. 3 is a block diagram showing an exemplary configuration of a randomnumber generation system according to a first exemplary embodiment ofthe present invention. As shown in FIG. 3, a random number generationsystem 100 of the present exemplary embodiment includes a generationmethod selection device 110, a rejection sampling device 120, and anaccumulation method sampling device 130.

FIG. 4 is a block diagram showing an exemplary configuration of thegeneration method selection device 110 according to the first exemplaryembodiment. As shown in FIG. 4, the generation method selection device110 of the present exemplary embodiment includes an interval uniformrandom number generation means 111 and a generation method selectionmeans 112.

The interval uniform random number generation means 111 has a functionof generating a uniform random number r in the interval [0, P].Furthermore, the generation method selection means 112 has a function ofselecting a random number generation method after comparing thegenerated uniform random number r and the probability p_({1-dim}).

In a case where the uniform random number r is probability p_({1-dim})or less, the generation method selection means 112 selects theaccumulation method as the random number generation method. After theselection, the generation method selection means 112 instructs theaccumulation method sampling device 130 to generate a random number.

In addition, in a case where the uniform random number r is larger thanthe probability p_({1-dim}), the generation method selection means 112selects the rejection sampling method as the random number generationmethod. After the selection, the generation method selection means 112instructs the rejection sampling device 120 to generate a random number.

FIG. 5 is a block diagram showing an exemplary configuration of therejection sampling device 120 according to the first exemplaryembodiment. As shown in FIG. 5, the rejection sampling device 120 of thepresent exemplary embodiment includes a uniform random number generationmeans 121 and a rejection determination means 122.

The function of the uniform random number generation means 121 issimilar to the function of the uniform random number generator 922. Thatis, the uniform random number generation means 121 generates a uniformrandom number u₁∈Z within a range of {k∈Z|−(t·s−a)≤k≤(t·s−a)}.

Next, the uniform random number generation means 121 generates areal-valued random number u₂∈R within a range of [0, φ(−a)]. The uniformrandom number generation means 121 inputs the generated u₁∈Z and u₂∈R tothe rejection determination means 122.

The function of rejection determination means 122 is similar to thefunction of rejection determiner 921 and the function of output device923. That is, rejection determination means 122 compares φ(u₁−a) andu₂∈R. When u₂≤φ(u₁−a) is satisfied as a result of the comparison, therejection determination means 122 outputs u₁∈Z as a random numberaccording to one-dimensional discrete Gaussian distribution the centerof which is not the origin.

When u₂>φ (u₁−a) is satisfied as a result of the comparison, therejection determination means 122 returns to the first step and executesthe same operation again.

FIG. 6 is a block diagram showing an exemplary configuration of theaccumulation method sampling device 130 according to the first exemplaryembodiment. As shown in FIG. 6, the accumulation method sampling device130 of the present exemplary embodiment includes a uniform random numbergeneration means 131, a search means 132, a storage means 133, and anoutput means 134.

The function of the uniform random number generation means 131 issimilar to the function of the uniform random number generator 913. Thatis, the uniform random number generation means 131 outputs a real valuex∈[0,1].

Furthermore, the function of the search means 132 is similar to thefunction of the searcher 912. That is, the search means 132 performsbinary search for z∈Z that satisfies φ(z−a−1)≤x<φ(z−a) from among thevalues stored in the storage means 133.

Furthermore, the function of the storage means 133 is similar to thefunction of the storage device 911. That is, the storage means 133stores the values of φ(0−a)/2, φ(1−a), φ(2−a), . . . , φ(t·s−a).

The output means 134 has a function of outputting the random numbersearched by the search means 132 as a random number according toone-dimensional discrete Gaussian distribution the center of which isnot the origin.

The random number generation system 100 of the present exemplaryembodiment is capable of generating a random number according to thediscrete Gaussian distribution on the multidimensional lattice at amemory cost for generating a random number according to one-dimensionaldiscrete Gaussian distribution by the accumulation method, and at acomputation cost close to the computation cost when each of the samplingprocesses is executed by the accumulation method.

Description of Operation

Hereinafter, operation of the random number generation system 100according to the present exemplary embodiment to generate a randomnumber according to one-dimensional discrete Gaussian distribution thecenter of which is not the origin will be described with reference toFIG. 7. FIG. 7 is a flowchart showing an operation of a random numbergeneration process performed by the random number generation system 100of the first exemplary embodiment.

First, the interval uniform random number generation means 111 generatesa uniform random number r in the interval [0, P] (step S101). Theinterval uniform random number generation means 111 inputs the generateduniform random number r to the generation method selection means 112.

Next, the generation method selection means 112 determines whether theinput uniform random number r is larger than the probability p_({1-dim})(step S102).

In a case where the uniform random number r is larger than theprobability p_({1-dim)} (True in step S102), the generation methodselection means 112 selects the rejection sampling method as a randomnumber sampling processing method (step S103). Next, the generationmethod selection means 112 instructs the rejection sampling device 120to generate a random number.

After instructed to generate a random number, the rejection samplingdevice 120 generates a random number according to one-dimensionaldiscrete Gaussian distribution the center of which is not the origin bythe rejection sampling method (step S104). After the generation, therejection sampling device 120 outputs the generated random number (stepS105). After the output, the random number generation system 100finishes the random number generation process.

In a case where the uniform random number r is the probabilityp_({1-dim}) or less (False in step S102), the generation methodselection means 112 selects the accumulation method as a random numbersampling processing method (step S106). Next, the generation methodselection means 112 instructs the accumulation method sampling device130 to generate a random number.

After instructed to generate a random number, the accumulation methodsampling device 130 generates a random number according toone-dimensional discrete Gaussian distribution the center of which isnot the origin by the accumulation method (step S107). The accumulationmethod sampling device 130 generates a random number by the accumulationmethod that uses a functional value defining one-dimensional discreteGaussian distribution the center of which is the origin.

After the generation, the accumulation method sampling device 130outputs the generated random number (step S108). After the output, therandom number generation system 100 finishes the random numbergeneration process. In a case where a random number according to thediscrete Gaussian distribution on the n-dimensional lattice isgenerated, the random number generation system 100 executes the randomnumber generation process shown in FIG. 7 n times.

Description of Effects

Since the value of the probability p_({1-dim}) is sufficiently largerthan the value of the probability p_({rjc}), the generation methodselection means 112 selects, in most cases, the accumulation method thatuses a functional value defining one-dimensional discrete Gaussiandistribution the center of which is the origin, as the random numbergeneration method. That is, when the random number generation system 100of the present exemplary embodiment is used, the sampling process thatis executed every time a random number according to discrete Gaussiandistribution on a multidimensional lattice is generated, is performed,with high probability, as the sampling process using the accumulationmethod.

Therefore, the memory cost required for the entire random numbergeneration process by the random number generation system 100 will bereduced to a value of a degree substantially the same as the memory costrequired when the random number according to the one-dimensionaldiscrete Gaussian distribution the center of which is the origin isgenerated by the accumulation method.

Furthermore, the computation cost required for the entire random numbergeneration process by the random number generation system 100 is reducedto a value close to the computation cost required when the samplingprocess executed every time is the sampling process performed by theaccumulation method.

Note that the random number generation system 100 of the presentexemplary embodiment may be implemented by a processor such as a centralprocessing unit (CPU) or a data processing device that executesprocessing according to a program stored in a non-transitory storagemedium, for example. That is, the interval uniform random numbergeneration means 111, the generation method selection means 112, theuniform random number generation means 121, the rejection determinationmeans 122, the uniform random number generation means 131, the searchmeans 132, and the output means 134 may be implemented by a CPU thatexecutes processes in accordance with a program control, for example.

The storage means 133 may be implemented by random access memory (RAM),for example.

In addition, individual components in the random number generationsystem 100 according to the present exemplary embodiment may beimplemented by a hardware circuit. As an example, the interval uniformrandom number generation means 111, the generation method selectionmeans 112, the uniform random number generation means 121, the rejectiondetermination means 122, the uniform random number generation means 131,the search means 132, the storage means 133, and the output means 134may be implemented by individual large scale integration (LSI) devices.Alternatively, they may be implemented by one LSI device.

Next, a summary of the present invention will be described. FIG. 8 is ablock diagram showing a summary of a random number generation systemaccording to the present invention. The random number generation system10 according to the present invention is a random number generationsystem that generates a random number according to a first discretedistribution, which is a discrete Gaussian distribution on aone-dimensional lattice the center of which is a positive value, therandom number generation system 10 including: a generation means 11 (forexample, the interval uniform random number generation means 111) whichgenerates a uniform random number between 0 and a first probability,which is a probability of a stochastic variable becoming a value withina predetermined interval in a positive range in the first discretedistribution; and a selection means 12 (for example, the generationmethod selection means 112), when a uniform random number less than orequal to a second probability is generated, the second probability beinga probability of the stochastic variable becoming a value within apredetermined interval in a second discrete distribution, which is adiscrete Gaussian distribution on a one-dimensional lattice the centerof which is the origin, selects, as a random number generation method,an accumulation method in which a functional value defining the seconddiscrete distribution is used, when a uniform random number greater thanthe second probability is generated, selects a rejection sampling methodas the random number generation method.

With such a configuration, the random number generation system canfurther reduce the memory cost and computation cost for generating arandom number according to discrete Gaussian distribution on amultidimensional lattice.

Furthermore, the random number generation system 10 may include anaccumulation method generation means (for example, the accumulationmethod sampling device 130) which generates a random number according tothe discrete Gaussian distribution on the one-dimensional lattice by theaccumulation method that uses the functional value defining the seconddiscrete distribution, the selection means 12 may instruct theaccumulation method generation means to generate a random number afterthe accumulation method is selected, and the accumulation methodgeneration means may generate a random number in response to theinstruction. Furthermore, the accumulation method generation means mayinclude a storage means (for example, the storage means 133) whichstores the functional value defining the second discrete distribution.

With such a configuration, the random number generation system cangenerate, by the accumulation method, a random number according todiscrete Gaussian distribution on a one-dimensional lattice the centerof which is not the origin.

Furthermore, the random number generation system 10 may include arejection sampling method generation means (for example, the rejectionsampling device 120) which generates a random number according to thediscrete Gaussian distribution on the one-dimensional lattice by therejection sampling method, the selection means 12 may instruct therejection sampling method generation means to generate a random numberafter the rejection sampling method is selected, and the rejectionsampling method generation means may generate a random number inresponse to the instruction.

With such a configuration, the random number generation system cangenerate, by the rejection sampling method, a random number according todiscrete Gaussian distribution on a one-dimensional lattice the centerof which is not the origin.

While the invention of the present application has been described withreference to the exemplary embodiments and examples, the invention ofthe present application is not limited to the above exemplaryembodiments and examples. Configuration and details of the invention ofthe present application can be modified in various mannersunderstandable for those skilled in the art within the scope of theinvention of the present application.

The above exemplary embodiments may also be partially or entirelydescribed as the following appendices, although this is not alimitation.

(Supplementary Note 1)

A random number generation system that generates a random numberaccording to a third discrete distribution, which is a discrete Gaussiandistribution on a one-dimensional lattice the center of which is anegative value, the random number generation system including: ageneration means which generates a uniform random number between 0 and athird probability, which is a probability of a stochastic variablebecoming a value within a predetermined interval in a negative range inthe third discrete distribution; and a selection means, when a uniformrandom number less than or equal to a fourth probability is generated,the fourth probability being a probability of the stochastic variablebecoming a value within a predetermined interval in a fourth discretedistribution, which is a discrete Gaussian distribution on aone-dimensional lattice the center of which is the origin, selects, as arandom number generation method, an accumulation method in which afunctional value defining the fourth discrete distribution is used, whena uniform random number greater than the fourth probability isgenerated, selects a rejection sampling method as the random numbergeneration method.

(Supplementary Note 2)

The random number generation system according to Supplementary note 1,further including an accumulation method generation means whichgenerates a random number according to the discrete Gaussiandistribution on the one-dimensional lattice by the accumulation methodthat uses the functional value defining the fourth discretedistribution, in which the selection means instructs the accumulationmethod generation means to generate a random number after theaccumulation method is selected, and the accumulation method generationmeans generates a random number in response to the instruction.

(Supplementary Note 3)

The random number generation system according to Supplementary note 2,in which the accumulation method generation means includes a storagemeans which stores the functional value defining the fourth discretedistribution.

(Supplementary Note 4)

The random number generation system according to any one ofSupplementary notes 1 to 3, further including a rejection samplingmethod generation means which generates a random number according to thediscrete Gaussian distribution on the one-dimensional lattice by therejection sampling method, in which the selection means instructs therejection sampling method generation means to generate a random numberafter the rejection sampling method is selected, and the rejectionsampling method generation means generates a random number in responseto the instruction.

(Supplementary Note 5)

A random number generation method to be executed in a random numbergeneration system that generates a random number according to a thirddiscrete distribution, which is a discrete Gaussian distribution on aone-dimensional lattice the center of which is a negative value, therandom number generation method including: generating a uniform randomnumber between 0 and a third probability, which is a probability of astochastic variable becoming a value within a predetermined interval ina negative range in the third discrete distribution; when a uniformrandom number less than or equal to a fourth probability is generated,the fourth probability being a probability of the stochastic variablebecoming a value within a predetermined interval in a fourth discretedistribution, which is a discrete Gaussian distribution on aone-dimensional lattice the center of which is the origin, selecting, asa random number generation method, an accumulation method in which afunctional value defining the fourth discrete distribution is used; andwhen a uniform random number greater than the fourth probability isgenerated, selecting a rejection sampling method as the random numbergeneration method.

(Supplementary Note 6)

A random number generation program executed on a computer that generatesa random number according to a third discrete distribution, which is adiscrete Gaussian distribution on a one-dimensional lattice the centerof which is a negative value, the random number generation programcausing the computer to execute: a generation process that generates auniform random number between 0 and a third probability, which is aprobability of a stochastic variable becoming a value within apredetermined interval in a negative range in the third discretedistribution; a third selection process that, when a uniform randomnumber less than or equal to a fourth probability is generated, thefourth probability being a probability of the stochastic variablebecoming a value within a predetermined interval in a fourth discretedistribution, which is a discrete Gaussian distribution on aone-dimensional lattice the center of which is the origin, selects, as arandom number generation method, an accumulation method in which afunctional value defining the fourth discrete distribution is used; anda fourth selection process that, when a uniform random number greaterthan the fourth probability is generated, selects a rejection samplingmethod as the random number generation method.

REFERENCE SIGNS LIST

-   010, 100, 910, 920 Random number generation system-   11 Generation means-   12 Selection means-   110 Generation method selection device-   111 Interval uniform random number generation means-   112 Generation method selection means-   120 Rejection sampling device-   121, 131 Uniform random number generation means-   122 Rejection determination means-   130 Accumulation method sampling device-   132 Search means-   133 Storage means-   134 Output means-   911 Storage device-   912 Searcher-   913, 922 Uniform random number generator-   921 Rejection determiner-   923 Output device

What is claimed is:
 1. A random number generation method selectingsystem that generates a random number according to a first discretedistribution, which is a discrete Gaussian distribution on aone-dimensional lattice the center of which is a positive value, therandom number generation method selecting system comprising: ageneration unit, implemented by a hardware including one or moreprocessors, which generates a uniform random number between 0 and afirst probability, which is a probability of a stochastic variablebecoming a value within a predetermined interval in a positive range inthe first discrete distribution; and a selection unit, implemented bythe hardware, when a uniform random number less than or equal to asecond probability is generated, the second probability being aprobability of the stochastic variable becoming a value within apredetermined interval in a second discrete distribution, which is adiscrete Gaussian distribution on a one-dimensional lattice the centerof which is the origin, selects, as a random number generation method,an accumulation method in which a functional value defining the seconddiscrete distribution is used, when a uniform random number greater thanthe second probability is generated, selects a rejection sampling methodas the random number generation method.
 2. The random number generationmethod selecting system according to claim 1, further comprising anaccumulation method generation unit, implemented by the hardware, whichgenerates a random number according to the discrete Gaussiandistribution on the one-dimensional lattice by the accumulation methodthat uses the functional value defining the second discretedistribution, wherein the selection unit instructs the accumulationmethod generation unit to generate a random number after theaccumulation method is selected, and the accumulation method generationunit generates a random number in response to the instruction.
 3. Therandom number generation method selecting system according to claim 2,wherein the accumulation method generation unit includes a storage unit,implemented by the hardware, which stores the functional value definingthe second discrete distribution.
 4. The random number generation methodselecting system according to claim 1, further comprising a rejectionsampling method generation unit, implemented by the hardware, whichgenerates a random number according to the discrete Gaussiandistribution on the one-dimensional lattice by the rejection samplingmethod, wherein the selection unit instructs the rejection samplingmethod generation unit to generate a random number after the rejectionsampling method is selected, and the rejection sampling methodgeneration unit generates a random number in response to theinstruction.
 5. The random number generation method selecting systemaccording to claim 2, further comprising a rejection sampling methodgeneration unit, implemented by the hardware, which generates a randomnumber according to the discrete Gaussian distribution on theone-dimensional lattice by the rejection sampling method, wherein theselection unit instructs the rejection sampling method generation unitto generate a random number after the rejection sampling method isselected, and the rejection sampling method generation unit generates arandom number in response to the instruction.
 6. The random numbergeneration method selecting system according to claim 3, furthercomprising a rejection sampling method generation unit, implemented bythe hardware, which generates a random number according to the discreteGaussian distribution on the one-dimensional lattice by the rejectionsampling method, wherein the selection unit instructs the rejectionsampling method generation unit to generate a random number after therejection sampling method is selected, and the rejection sampling methodgeneration unit generates a random number in response to theinstruction.
 7. A computer-implemented random number generation methodselecting method to be executed in a random number generation methodselecting system that generates a random number according to a firstdiscrete distribution, which is a discrete Gaussian distribution on aone-dimensional lattice the center of which is a positive value, therandom number generation method selecting method comprising: generatinga uniform random number between 0 and a first probability, which is aprobability of a stochastic variable becoming a value within apredetermined interval in a positive range in the first discretedistribution; when a uniform random number less than or equal to asecond probability is generated, the second probability being aprobability of the stochastic variable becoming a value within apredetermined interval in a second discrete distribution, which is adiscrete Gaussian distribution on a one-dimensional lattice the centerof which is the origin, selecting, as a random number generation method,an accumulation method in which a functional value defining the seconddiscrete distribution is used; and when a uniform random number greaterthan the second probability is generated, selecting a rejection samplingmethod as the random number generation method.
 8. A non-transitorycomputer-readable capturing medium having captured therein a randomnumber generation method selecting program executed on a computer thatgenerates a random number according to a first discrete distribution,which is a discrete Gaussian distribution on a one-dimensional latticethe center of which is a positive value, the random number generationmethod selecting program causing the computer to execute: a generationprocess that generates a uniform random number between 0 and a firstprobability, which is a probability of a stochastic variable becoming avalue within a predetermined interval in a positive range in the firstdiscrete distribution; a first selection process that, when a uniformrandom number less than or equal to a second probability is generated,the second probability being a probability of the stochastic variablebecoming a value within a predetermined interval in a second discretedistribution, which is a discrete Gaussian distribution on aone-dimensional lattice the center of which is the origin, selects, as arandom number generation method, an accumulation method in which afunctional value defining the second discrete distribution is used; anda second selection process that, when a uniform random number greaterthan the second probability is generated, selects a rejection samplingmethod as the random number generation method.